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ABS'TRACT 


Although parts of this paper may be applicable to various 
other location finding systems, radio direction finding was 
selected for development because of the extensive literature 
available on the subject. Using a modification of the navi- 
gation coordinate system, equations in spherical trigonometry 
are developed by analogy with graphical solution procedures. 
Ellipticity of the earth is corrected for by use of geocen- 
tric latitude in all calculations. The least squares solu- 
tion is found by a combination of second-order gradient 
searches and Golden Section Searches on the corrected sphere. 
A Chi-Square contour confidence region is presented with 
defense against the empty region arguments and with verifi- 
cation of its reliability of target containment. Sensitivity 
of location error to network layout, number of measurements 
available and intercept geometry factors is examined by 
computer simulation, demonstrating the predictive value of 
these statistics. It is shown that the correlation between 
location error and the size of a reliable confidence region 


is too minimal to have predictive value. 
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INTRODUCTION 


Various procedures have been described for estimating an 
emitter location from a set of measured signal azimuths 
reported to a control station by dispersed direction finding 
(DF) sites. Such procedures range upward in complexity from 
visual estimates off a plot of great circle arcs on a 
gnomonic projection of the operating area of potential 
targets. Other procedures exist for describing confidence 
criteria about the location estimate. The present paper was 
motivated by an interest in confidence regions and in loca- 
tion error analysis. Such studies require either that the 
procedure used be free of approximating assumptions, or that 
the variability of the approximations be included in the 
confidence criterion. The former course was preferred for 
simplicity and consistency of estimates. This paper presents 
an emitter location estimate procedure based only upon known 
DF site coordinates and signal azimuths. The procedure 
described in sections II and III converges quickly to a 
least squares solution and does not require an initialization 
point estimate of the location from external sources (Pope, 
1971). The procedure does not require a rejection of out- 
ters (McCalla, 1970) in order to converge to the least 
Squares solution. The approach is free from assumptions in 
the categories that the earth is flat in the solution region, 


that the earth is a perfect sphere, that great circle arcs 





are well enough approximated by straight lines of position, 
and that azimuth errors are well enough approximated by 
parallel displacements of such lines of position (Daniels, 
ESSI). 

Section IV describes a Chi-Square contour confidence 
region procedure. Factors from intercept geometry are 
discussed in section V. Section VI is a description of 
procedures used to simulate the problem for the purpose of 
the location error analyses described in sections VII and 
VIII. An alternative to confidence regions for predicting 


lion accuracy 1s discussed in the concluding section. 





IN  IIUNPICATSDESCREPTION OF SPATIAL RELATIONSHIPS 


Locations and spatial relationships of interest in direc- 
tion finding and target fixing Can be expressed completely 
in spherical trigonometry. The advantage of such expression 
is that each mathematical step has its simple graphical 
counterpart. This paper contains notational and descriptive 
departures from other direction finding models based on spher- 
ical trigonometry. Such notational changes follow those used 
in geodesy for describing locations and spatial relationships 
and those used in textbooks on spherical trigonometry for 
developing equations to express such relationships.  Fquations 
are restricted to sine and cosine functions only in order to 
simplify calculating first and second order derivatives of 


somewhat lengthy equations. 


E GLOBAL COORDINATE SYSTEM 

The coordinate system is a slight modification of the 
System used in E Longitude, denoted À, is measured 
from Zero at the prime meridian eastward through 360 degrees. 
Marrtude, denoted ¢, iS measured from zero at the equator 
northward through ninety degrees and southward through minus 
ninety degrees. Azimuth, denoted 6, is measured from zero 
at true north clockwise through 360 degrees. Such a coor- 
dinate system is sufficient to represent any point on the 
globe as well as the great circles which represent ideal wave 


propagation paths. Mathematical descriptions of spatial 





relationships must account for the longitude discontinuity 


at the prime meridian and for the azimuth discontinuity at 


fue north. A great circle arc with azimuth angle indicated 


is sketched on the coordinate system in Figure 1. 





Global Coordinate System 


FIGURE 1 


Bee GEOCENTRIC LATITUDE 


Geographical latitude, used in navigation and mapping, 


is based upon astronomical observations and mathematical 


Hedaetxons that do not correct for the ellipticity of the 


earth. 
greater 
equator 
earth. 
applied 


version 


Consequently, geographical latitude is slightly 

in magnitude than the arc on a meridian from the 

to a location with vertex angle at the center of the 
Rules of spherical trigonometry can be more accurately 
to describe spatial relationships following the con- 


from geographical latitude to geocentric, or earth- 


10 





centered, latitude. Assuming that waves propagate along 
great circle arcs, then no information is lost since such 
paths are geometrically independent of radii to the earth 
enter from both the signal emitter and from the DF site. 
The maximum discrepancy between the two types of latitude 
occurs at +45 degrees latitude and is about equal to 11.5 
minutes of arc. Although the conversion equation appears 
in books on geodesy, a derivation could not be found, so the 
equation was derived for reference purposes. 

In Figure 2, a and b represent the semi-major and semi- 
minor axes of the ellipsoid. The point P (xi, yi) represents 
any point on the ellipsoid. The ine OP connects P to the 


origin at O, the earth center. SP is perpendicular to the 


D 






HaUalioll Oi ellipse: 


2 2 


E 
2 
a 





N S a 


Relationship Between Geographical and Geocentric Latitude 


FIGURE 2 
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equatorial plane, the x-axis. NP is the normal line to the 
ellipsoid at P, and the angle PNS is equivalent to geographi- 
ЕСИ Таб Е1ае. Angle PON is equivalent to geocentric latitude, 
denoted Y. Ellipticity of the quadrant is exaggerated in 
Order to distinguish the two angles. 


From Figure 2, it can be seen that 
СГ ООС tan WU, ana that SP = NS tand, 
so that 


NS 


tany = OS tano ZE 


From the equation of the ellipse, 
(a^ - xD) (2-2) 


The equation of the normal line NP is given by 


1 


Y-Yi^ - fry 


(IFE X) E) 


Differentiating equation (2-2) to solve f'(x)), 


eni 
С 


' NEM S 2 
i (х) E a 


2 
E (a - ху) 


Substituting in equation (2-3) and rearranging terms, 


1 
2 


E 2.2 E 
2 
a yi 
= — (x - xj), from equation (2-2). 
bx 


Then at y = 0, y, can be cancelled from both sides, and 


®2 





e 
Xi a2 
From Figure 2, ху T ас у = O, x = ON, so that 
Xi — x = NS 
NS _ b“ 
OS 72 
a 


Substituting this last result in equation (2-1), 


2 


&апф = 25 tan (2-4) 
a 


The conversion is usually stated as a function of earth 


eccentricity, 


es => le =b) 
Then 
Ee B a* - s 
Е 2 
а 
IE: 
2 
a. 
Therefore, 
Pec 2 | 
E l- e, and it follows that 
a 
ú 2 
tany = (1 - e^) tanó (255) 


Using e = 0.081813 from the International Fllipsoid of 


Reference (reference 3), 


tanyp = 0.99327733 tanó (226) 





C. THE POLAR TRIANGLE 

The mathematical description of spatial relationships is 
based upon a polar triangle constructed from the great circle 
arcs (measured by the angle subtended by the arc at the 
earth's center) connecting a DF site, any other point and the 
north pole. Such a triangle is sketched in Figure 3. The DF 
is located at vertex B. Vertex C is at any other point. 
Vertex A is at the north pole. The arc c is the complement 
of the geocentric latitude of the DF site. Similarly, the 
arc b is the complement of the geocentric latitude of the 
other point. Angle A is equal to the difference in longitude 
between the DF site and the point. The dotted arc included 
on Figure 3 corresponds to a signal azimuth with angle indi- 


cated by 80. 





Polar Triangle with Signal Azimuth 


FIGURE 3 
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Angle A corresponds to the smaller of two possible polar 
triangles connecting any two points and the north pole. 
Usually, there is no problem in deciding which polar tri- 
angle to choose. However, when this triangle contains the 
prime Meridian, the discontinuity at longitude 360 = longi- 
tude 0 results in assigning the value of the vertex of the 
larger polar triangle when the longitude difference is taken. 
The mathematical model can be forced to correspond to the 
spatial relationships by adding the proviso that if the 
difference in degrees of longitude is greater than 180 
degrees, then let A equal 360 degrees minus the difference. 
This restricts the size of angle A to one hemisphere. The 
absolute value of the difference can be taken in case the DF 
site is east of (greater longitude than) the point. These 
relationships are summarized below with subscript s desig- 


nating the DF site and subscript t designating any test 


Point. 
ce = 90 - y. (2-7) 
BE 90 u (2-8) 
Ma Es 180 
pu (2-9) 
360 - I^. - Al, a E JA, ~ A | > 180 


mE orne Azimuth Solution 
The value of arc a connecting the DF site and test 


pountorsoweiven by the law of cosines for sides, 


COS a = cos b cos c * sin b sin c cos A (2-10) 
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There is no quadrant ambiguity since the cosine takes nega- 
tive values in the second quadrant and arc a is always less 
than or equal to 180 degrees because of the restriction on 
angle A to one hemisphere. Since the cosine of angle A is 
EwHnterest in determining arc a, both the absolute value 
sign and the condition can be removed by applying a reduction 


formula. From equation (2-9), 


cos |A, - X | 
cos A = : 
E EE (A, = A] > 180 
Eu through reduction, 
cos (360 - A, - UD = cos JA, - Х| 
The absolute value sign can be dropped since cos (A, = A) = 
cos (N. - AL). From this result and equation (2-10) 
Cos a= cos b cos ec rsin b sin c cos (A, - е, (2-11) 


The vertex angle B at the DF site can be found by applying 


Ane Of the laws of cosines for sides, 


cos b —- cos c cos a 


Sin C sina (2-12) 


cos B = 


Equation (2-12) is valid for test points either east or west 
of the DF site at vertex B. 

In order to make comparisons with a signal azimuth, 
u must be converted to an azimuth from true north, 
denoted B. Visually, the conversion is based on whether the 
test point is eastward or westward of the DF site with res- 


pect to the proper polar triangle.  Mathematically, the 


l6 





conversion is based on the magnitude of the longitudes of 
NNF ete ana the test point, but must also account for 
the discontinuity at the prime meridian. There are four 
general cases with two distinct outcomes as illustrated in 
Figure 4. 

2. Point Azimuth and Signal Azimuth Difference 

From a given DF site, the difference between the 

azimuth of any test point and a signal azimuth is the smaller 
Of the two possible outcomes, depending upon whether the 
difference is measured clockwise or counter-clockwise. The 
sun h difference is denoted o, and the azimuth discontinuity 
at true north requires that the computation be conditional. 
With a clockwise rotation for the azimuths, the smaller 
angle of azimuth difference will be positive if 6 is 
leading and negative if 8 trails 0. The possible outcomes 


Enc their conditions are as follows. 


a. If |B - 6| < 180, then w = 8 - Ө 

b. If [8 - 6| > 180 and 8 > 6, then w = B - 0 — 360 

c. If |8 - 0| » 180 and 8 « 6, then w = 8 - 6 + 360 
| (2-14) 


The sets of conditions on computing B (B) in equation 
(2-13) and in computing u(8,0) in equation (2-14) yield six 
possible outcomes for computing azimuth difference as a 


function of the vertex angle B and signal azimuth 6. 
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ge sand Ay = N < 180, then B = B 





< 180, then B = 360 - B 


б МЕЛ ЖА апа hz EE 


ë Ar > 180, then B = B 


ie 


д 


c В < À and hy - ic > 180, then B = 360 —- B 
Azimuth and Vertex Angle 


FIGURE 4 
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BE > 360 = 9 


БЕ СО O 


u = 
со = о 
XE O 
ZO = в 


Two utility variables may be introduced in order to express 
azimuth difference in a constant form with values dependent 


upon the existing conditions. 
Ww = pB +q - 8 (2 15) 


ШЫ ра1г of utility variables can take on any one of the 


following sets of values 


(150) 

(1,-360) 

(360) 
(p,q) = 

17360) 

(-1,0) 


=, 720) 


Concluding development, angle B can be expressed as a differ- 
entiable function of sines and cosines of known values of 
longitude and geocentric latitude. From equations (2290) 

and (2-12) and the application of an angle-difference rela- 


Elon, 
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p = cos +} (co b=) Cos с [cos mesos tes + Sin b sin cC 


(cos A, cos A + sin A, sin JJ / C 


А ~] : : : 
sin (cos [cos pecos OC t sin b sin c (cos A, cos 


+ sin A, sin xl 


Substituting the values for arcs b and c from equations (2-7) 
and (2-8) and using the reduction formulas cos (90 - x) - 


sin x and sin (90- x) = cos x, 
_ -1 |f... 65 À . 
p COS (6 UM sm Y. [sin V, sin V. + cos Y, 


COS Ve (cos Ar cos À + sin AL sin 1911) / [cos Yo 


: -1 3 В 
Sin (COS n V. sin Yo + cos Y, cos Ve (cos A. 


cos A. + sin A, sin a] 


By multiplying the term sin Vo through the interior brackets 
of the numerator, writing the first two resulting terms as 
sin Vi (1 - sin^y ), which is equal to sin Y, səs 


cos Yo can be cancelled from both numerator and denominator. 


_ = а 
В = cos [е Y, COS Ú. sin V. COS pe (cos AL 


: ; ; -] Р 
COS À + sin AL sin | /| sin (cos | sin Y, 


sin V. + cos Ú+ cos V. (cos AL COS À 


+ sin A, sin ))]) | (2-16) 
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IAE TEA SOLUTION 


In Figure 5, each of four DF sites have measured azimuths 
of a signal from the same emitter at an unknown location. 
The signal azimuths are plotted for illustration as dotted 
ARES. A test point, PO, V), is established, and azimuths 
Beem ach DF site to P are drawn in solid arcs. A better 
fit can be found by moving P somewhere within the region. 
bounded by the dotted arcs. 

Given a number n of DF sites, each reporting the azimuth 


of a received signal with respect to its own location to a 


net control, the fàx problem can be viewed as finding the 





Azimuth iar ferences 


FIGURE 5 
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point on the globe which minimizes the sum of some function 


of the azimuth differences We TOT 1 l210. 


OBJECTIVE FUNCTION 

Since the azimuth differences are signed, minimizing 
their sum may produce a solution that lies outside the region 
bounded by arcs corresponding to signal azimuths when large 
positive and large negative azimuth differences cancel each 
other out.  Minimizing the sum of the absolute values of 
azimuth differences avoids this problem of positive and 
negative differences; however, since the partial derivatives 
of absolute value are not computationally tractable, solution 
techniques are limited. Minimizing the sum of the squared 
azimuth differences provides an objective function that is 
differentiable and reflects absolute magnitude. Therefore, 
the solution of the fix problem is taken to be that pair of 


Coordinates Vir V.) which will 


I a ea 


Minimize 


Pile tees 


D 


B. SECOND-ORDER GRADIENT SEARCH 

Dropping Ene t subscript to simplify exposition, the 
objective function, equation (3-1), can be represented as 
f(A,U). Given an initial test point (A,W), a truncated 
May lors series for a function of two variables can be used 
to approximate the functional values of other test points 


fet AA, Y + Ay). 
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f(A + AA, w+ Ab) 2 £(A,U) + Ay 2 N ED 


IE (1,4) 


2 
DES E (AD) 1 
2 343g ^ 20V) 


2 3“f (X ,1p) 
3À 2 


д 


T = (AX) + AMY 


Expressing the approximation in matrix form and letting 


Enperscript t indicate the transpose, 


ЭЕ (Ауф) 
A + Ad Е JA 
f - £(A,U)* + (АА, ду) 
Y + Ay IF (ALY) 
3j 
E (AY)  23^£ (Ou) 
7 p Seh! AA 
+ 5 (AA, Ay) 
ea) E,W) au 
AD e 


Letting Р, denote an initial test point CPUS letting R 
denote the vector (ЛА ,ду) Ё, letting V denote the gradient 
Or £ (P) and letting y denote the Hessian matrix of f at 


Р., the approximation can be expressed as 


f(D. + R) = £(P.) « R've(P.) + | Rv f(E.) R 
j j j j 
Then, 
Е an E + > O EM 
£(P, * R) - £(P,) 3 RVE(B) + 7 RVOE(B,) R (3-2) 


Starting from an initial test point at Р, the search pro- 
cedure can be understood as a technique to find a vector R 
such that Ё (Р. + К) is the minimum of the objective function. 


Such a vector R also minimizes the left hand side of equation 


E25. Since the amada rata forma on the right hand side is an 


23 





approximation of the value of the left hand side, it must 

be near a minimum value when the left hand side is minimized. 
The approximation works both ways, so that a vector R which 
minimizes the right hand side must approximately minimize 

the left hand side. The right hand тае is in a form that 
ds to analysis of its extreme values with respect to R 
at a fixed Р.. One of the conditions necessary for the 
quadratic form to be at a minimum is that its gradient be 
equal to zero. Therefore, it is necessary to find a vector 
ПИЕСА fulfills this condition. Taking the gradient of 


Esustion (3-2) with respect to R, 


ее ык z 
ZI VER) + HR VER.) 1 
— v. [stes 115622,5. 5 
= А; А) + = куе (Б.) R 
- 1 22.5.5 1562,5 
ل‎ Ето. 
( 3! > ( 3) > ( 3) 


_ Dee - 
VEE FP) R 
3) J 
Mina the gradient of the quadratic form at R equal to zero, 
= 2: E = 
Mu P cv f(P-) R = 0 
J J 
Then, 
V^£(P.) R » -V£ (P) 
J J 
And, providing that ve (BS) 1s invertible, 


= Dn m = 
В = -—- V «y BI (3-3) 
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R specifies both a direction and magnitude in which to 
move toward an extreme value of the objective function. 
This consequence of the Taylor approximation can be improved 
by applying a scalar coefficient, denoted 6, to R in order 
to adjust the magnitude and, if R points to a maximum value, 


the direction. The next test point is then given by 


. = Бр + 6R (3-4) 


Somoining equations (3-3) and (3-4), 


ET [0220.3] КЕСА СЕ, 
The search procedure can be terminated when further 
improvements are of little consequence, most readily deter- 
mined by comparing the distance between locations on succes- 

Sive iterations. 


In order to solve equation (3-5), it is necessary to 


compute the inverse Hessian of the objective function. 


3^£(P) — 97£ (3) 


"e = key т 3y* элэ 
ү f (P) = Mo 2 а 2 e 
det|v^r(p)| | s^r c) ae 

эд?‏ ل3۸3 


where the determinant of the Hessian 


IN 3:sB 6: E 
det |v #(Р)] = E n E 
SA 3U дА дф 
Substituting the inverse Hessian in equation (3-5) and expand- 
ing the equation by matrix multiplication of the inverse 
Hessian times the gradient yields a new test point with ó the 


only unknown. 
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= о Ы 
3! 3f (P4) д £ (P4) df (Pj) 


P 
TE 3۸ JAI 90 
p 2. ,3 2. ,= 
ИО ЕУР ТИ OER) 


2 S 3۸234 


—- <  —— PP س‎ шш n 


9 À 


"Jl 
II 

iu) 
l 

О» 


Jt] j 

34£(P.) £ (Fj) 3^ f (P.) 0£ (F4) 
2 50 TOW 3À 
12 


9A 
2 2 - 2, 
у? Ё 9^ f (P; 
317 9р? 3۸34 


(9:59 


The solution of equation (3-6) in turn requires the cal- 
Station of all first and second partial derivatives of the 


objective function. 











JA 
n 9B; 
= 2 2p; (p;Bi + ei = ө.) БШ (3-7) 
а 
ОП 1аг1у, 
= Ä oB. 
Жр Л 7 i a 
2 2 
2 £z n oB. o B. 
9 f (P) 1 1. 
GONE - 5) 2p + 2p. (p.B. + q.- 0.) 
(3—9) 
2 2 
2 - ,= n oB. o D: 
9 f (P) 2 1 т 
اا‎ = E 2p? |——| + 2p. (p.B.+q. - 6.) 
9р? i=] 1 OU a 1 1 TE 
(3-10) 
pM n SE OE 32B. 
EP) O 2 2 m > n 
EE Pi sy ?Pi PP; * d; 91) 315g 
(3-11) 
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Solving equations (3-7) through (3-11) requires the 
Nena thy procedure of taking all first and second partial 
derivatives of equation (2-16), the expression for angle B. 
Solutions to these equations are contained in the FORTRAN 


Fix Program included in this paper. 


GOLDEN SECTION SEARCH 


From equations (2-15) and (3-1), the objective function 


can be written, 


(p,B, + a, - 8,)^, E) 


Tl 
. IL 


Minimize 
i=l 
the sum of the squared azimuth differences. 


Letting x denote the indicated move in longitude and y 


denote the indicated move in latitude given by equation (3-6), 


R = (AA, AV) 


(x,y) © 


N 


Then angle Bis given in equation (2-16), becomes at zem 


PE o БАЗЫ cosV; - sinV, rn 


i 
е е a -1 
À- ; -= ; 
[eos 3 ӧх) cosA, +sin(A, ôx) лаа ]) / sin [cos 


[sin (v,-5y) sinW, * cos (V; - Sy) cost), [cos (A, - 5x) 


| (3-13) 


where subscript j denotes test point coordinates and sub- 





cosh. + sin (A, - 6x) sind, | 


dript i denotes DF site coordinates. Since test point 


2] 





coordinates are arbitrarily known and DF site coordinates 

are constant for the problem, the substitution of equation 
EL) into equation (3-12), the objective function, results 
in the problem of minimizing a function with a single unknown, 


ШОО ле scalar coefficient. 
аы 2 
Minimize ОЕ x w; (ò) (3-14) 


The objective function could be evaluated at incremental 
K s Of ô, say over the interval (a,b). Then that 6 which 
yields the smallest functional value would be chosen as the 
scalar coefficient to adjust the move indicated by the 
second-order gradient NEC The Golden Section Search 
(reference 1) provides a fast and methodical approach to find 
this minimizing ó. Figure 6 illustrates the problem; however, 
the actual curve representing f(ó) is unknown at the outset 


of the problem. 


— 


Reducing the Search Interval 
FIGURE 6 
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Based on the assumption that f(6) has a local minimum on 
(a,b), the objective is to find that minimum without evalu- 
Bene the function at too many points in (a,b). The Golden 
Section Search systematically eliminates portions of the 
interval by comparing functional values of pairs of points 
EE 3,5). From Figure 6, £ (65) is less than £ (94); there- 
fore, the portion (a,64) can be eliminated inesearehing for 
the minimum, reducing the search interval to ($1,b). A new 
point then is selected for comparison with © (6), resulting 
in a second elimination of whichever end re el corresponds 
to the higher functional value. 


ЕЙ 


Letting I. represent the length of the n search inter- 


val, 


Two different values of 6 are tested at each trial. By 
letting the two values be symmetric on I, the length of 


ia is the same regardless of which end interval is dropped. 


Letting o denote the relative position of $4 in Ii, 
go = (6, - a)/(b - a) 
Then 


$1 = G(b - a) ta (3-15) 


Since 6, is symmetric with $i: the relative position of Š in 


I; is given by (1 - 0), so that 


1 - cç = ($5 EMI za 
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and 


$ JM) (ta) + a IL E 


Provided that $ yielded the minimum functional value, as in 
Figure 6, then $ becomes the first element of the symmetric 


pair in I, and its relative position in I, is given by 


2 


(Š 


O ASADAS a) 


2 


| 


2-0) - а) +а ~ 6/1 - о) (Б E 


from equation (3-16), so that 


6, = (1 - x S apama sls (3-17) 


From equations (3-15) and (3-17), it follows that 


HE uoo a) a 
с = (1 - o)? (3-18) 


Mon the other hand, $1 yields the minimum functional value, 
then $1 becomes the second element of the symmetric test pair 


in I5, and its relative position in I, is given by 


2 


1-o = (6, - a)J/(1 - 0)(b - a) 


1 


So that 


DO e m‏ - 1) = وة 


as in equation (3-17), and leading to the same result as in 
equation (3-18). 
Brom equation (3-18) 


Bey 


which can be solved by the quadratic formula, so that 
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a 
| 


[з + (9 - 4) *]/2 


07382 


Ц 


En rtuting this result in equations (3-15) and (3-16), 


Or 
| 


0:382 (b - a) +a (3-19) 


5 


0.618 (b - a) + a (3-20) 


Generalizing, let b. and a, represent the endpoints of the 


Ben Search interval. Then 
$1 DASS (b. _ an) + a, (3-21) 
$ = 0.618 (b. -- a,) + a, (3-22) 


/ 


l' cS б, Or $. 1s the 


survivor from the preceding iteration, so that a comparison 


After the first comparison on I 


can be made after evaluating only the new test point symme- 
СУУ ЕБ the survivor. The Golden Section Search terminates 


when b. - a < g, where є пау ре tailored to the magnitude of 


n 


mie indicated move. 


The value of ó which minimizes f(6) over (a, ,bi) is 


entered in equation (3-6), resulting in a new test point 


t 


5 E t Е 


D. FORTRAN FIX PROGRAM 


The equations presented in this paper were coded in 


FORTRAN. The resulting program has been successfully 
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EXccuted on an IBM 360/67 computer system in the W. R. Church 
Computer Center at the Naval Postgraduate School. 
Program Description 

nbecgo5ogramecontarns an identification number, longi- 
tude, latitude and name for each DF site in an imaginary 
network. Inputs to the program are DF site identification 
number and signal azimuth. After DF site names and coordi- 
ees mT Correlated with input site identification numbers, 
Nun trigonometric functions of site coordinates are 
computed and stored in an array. Searching for the minimum 
sum of azimuth differences begins from (0,0) with a Golden 
Section Search on longitude using (360,0) as the indicated 
move. The result is a minimizing longitude on the equator 
from which a Golden Section Search is conducted on geocentric 
Eide with indicated move (0,180). From this minimizing 
geocentric latitude, another Golden Section Search is con- 
ducted on longitude with (180,0) as the indicated move. 
Such alternating searches on longitude and geocentric lati- 
tude continue until diminishing returns are indicated by a 


new position being less than 500 nautical miles from the 


previous minimizing position. The initial interval for the 
Golden Section Search is set at (-1,1), the endpoints for 
equations (3-19) and (3-20). Subsequent endpoints are 


modified by the search results, producing values for equa- 
tions (3-21) and (3-22). 
The program then executes second-order gradient 


searches refined by Golden Section Searches on indicated 
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moves until the most recent move is less than some input 
termination value, say one-tenth of a nautical mile.  Indi- 
cated moves, given by equation (3-6), are unbounded and 
pent occasionally go off the coordinate system of Figure i. 
Scaling any such wild moves to the coordinate system by 
modulo arithmetic destroys both the direction and magnitude 
of the vector. The indicated moves are kept in bounds by 
conversion to a unit vector in radians until such moves are 
of magnitude less than 30 nautical miles, such a small step 
imercating that the solution is at hand. By converting the 
maampcated move to a unit vector prior to refinement by a 
Scalar coefficient from the Golden Section Search, moves 
from the second-order gradient search are limited to ahout 
3440 nautical miles per iteration. 

Search output on succesSive iterations is optional, 
depending on an operator entry. Unless search progress is 
requested, searching terminates after thirty iterations even 
if the latest step is greater than the input termination 
value. After every tenth second-order gradient search, if 
Such a. high number should be required, the search effort is 
jolted by conducting alternate Golden Section Searches on 
longitude with indicated move (90,0) and latitude with indi- 
cated move (0,90). Such Golden Section Search interruptions 
only yield a new test point when they reduce the current 
value of the objective function. 

То lO SOIU ION COUCPDUL, geocentric latitude is 


converted to geographical latitude using equation (2-6), and 
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the solution is given in degrees and minutes of longitude, 
east or west, and latitude, north or south. 
2. Results 

The program has worked successfully in over ten 
thousand runs against simulated targets spread randomly over 
the globe with various networks of DF sites and various 
numbers of sites intercepting a given target. In actual 
Æ etice, prior information about a target's area of opera- 
night produce a better starting point than (0,0); how- 
ever, the time required to type such starting coordinates 
on a keyboard for input would exceed the total time required 
for the solution starting from (0,0), which averages less 
than a second and includes some diagnostic output useful in 
analyzing results. Figure 7 illustrates a representative 
set of search iterations with Golden Section Searches 
plotted as dashed lines and second-order gradient searches 


Solid. 


(0,0) 


Searching 


FIGURE 7 
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BN ARGEBEDOCATPIONSCONFIDENCE REGION 


That location yielding the minimum sum of squared azimuth 
differences, or least squares solution, would coincide with 
the location of the signal emitter if all measured signal 
azimuths with respect to the DF sites were free from error. 
Since then all great circle arcs corresponding to signal 
azimuths would intersect at a point, the minimum value of 
the objective function would be equal to zero, as is always 
the case when only two signal azimuths are available to 
compute a fix. Because the sites in a DF network do not 
produce errorless measurements, solutions rarely coincide 
with the emitter location, and there is the question of how 


close the emitter might be to the least squares solution. 


E SIGNAL AZIMUTH ERROR DISTRIBUTIONS 

Supposing that signal azimuth measurements made at indi- 
vidual DF sites are not biased in any of the possible direc- 
tions, then it follows that signed errors have a mean value 
of zero. If such measurements were found to be biased in 
certain directions during routine calibration trials, then 
the amount of bias could be applied as a correction to 
subsequent measurements, producing an estimated mean value 
of zero.  Lacking any reason to assume a different distribu- 
tion and following the lead given in the literature of 


u Lion finding, the error distribution can be taken to 
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be normally distributed with mean zero and variance т, 
assumed known and independent from site to site. 

The error distribution can be made to be symmetric 
around zero by following the sign convertion that clockwise 
errors are negative and counter-clockwise errors are posi- 


tive in sign. 


B. M ROBJECTIVE FUNCTION VALUE AS A CHI-SQUARE STATISTIC 
From standard textbooks on probability theory (for 

example, reference 2), it can be verified that if a random 
variable is normally distributed with mean zero and variance 
E then the square of that random variable is governed by 
Eie Gamma distribution with parameters (1/2,1/20%). катет 
more, the sum of n independent Gamma distributed random 
variables divided by the variance of the underlying normal 


distributions is a Chi-Square statistic with n degrees of 


K Om. Hence, where w denotes the azimuth difference, 





2 
0;7x, N(0,0:) (4-1) 
е (4-2) 
26 
ak 
n 
y Zur, (e 
i-1 01 


Since the normal distribution here referred to applies 
to the azimuth differences between the measured signal 
azimuth and the azimuth of the unknown emitter location, the 
ш. refer to actual and unknown errors rather than calculated 


differences between solution azimuths and signal azimuths. 
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Equation (4-3) therefore specifies the distribution of 
values of the objective function (equation (3-12)) at unknown 
emitter locations under the assumption that azimuth errors 


are normally distributed with known variance. 


See CONFIDENCE LEVELS AND REGION BOUNDARIES 

Since the sum of squared azimuth errors divided by 
Variances is a Chi-Square statistic under the previous 
assumptions, there is a specific probability, denoted a, 
that this sum will be less than or equal to any given value. 
Benavalently, for a fixed probability a, or confidence 
level, there corresponds an upper bound, denoted doe 


which depends upon the degrees of freedom as well as the 


eonfidence level. 


а (4-4) 


ime Choice of confidence level, a, 1S arbitrary, The 
corresponding value for aa can be found in Chi-Square 
tables. 

For a given fix ler the value of the sum of least 
squares, denoted f*, is the minimum possible value of the 
Chi-Square statistic for the observed signal azimuths. 
Furthermore, by evaluating the objective function at suc- 
cessive points outward in any direction from the least 
squares solution, it can be seen that the farther away the 


point evaluated up to the antipode, the higher the objective 


Kul On value. The location corresponding to f* can 
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therefore be regarded as a point in a region which contains 
all the possible emitter locations for which the value of 
the objective function is less than or equal to an upper 


pound C given the set of observed signal azimuths. 


(n,a) ” 
Such a region is bounded by all locations for which the 
objective function has a value equal to OUT and it can 
be stated that, for the observed signal azimuths, there is a 
probability a that the bounded region will contain the 
unknown emitter location. 

Daniels and Vajda (1951) challenge the utility of such 
an approach on the grounds that at given confidence levels 
the corresponding region may not exist with positive 
probability if n is greater than two. Their observation is 
a consequence of the fact that except for n equal to two 
Gi the great circle arcs corresponding to signal azimuths 
seldom intersect at a single point thereby yielding f* 
equal to zero. 

Prior to the availability of the high-speed digital 
computer, the practical problems of producing such a region 
were certainly more striking than theoretical objections. 
Even so, the Daniels/Vajda criticism was more than casual, 
which may indicate a requirement for additional arguments 
to justify the use of such a confidence region. | 

Suppose that f* is equal to k for some least squares 


solution, where 


rg 
ПКБ 


E =, (4-5) 


E 
uu c? 
i 
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and the probability p may be determined by interpolating 

with the value k in a table of Chi-Squared values. Daniels 
and Vajda are concerned that if the confidence level a had 
been initially specified at some value less than p, then 
there exists no region of locations having objective function 
values less than or equal to k, k being the minimum possible 
value on the globe for the given signal azimuth observations. 
There are at least two responses to such an occurrence and 
the accompanying empty region objection: (1) Increase the 
confidence level whenever the region does not exist at the 


maior confidence level so that Ci is greater than k, and 


no) 
(2) Routinely specify a confidence level of sufficient mag- 
nitude to make the occurrence of an empty region highly 
improbable. 

Either of the above alternatives insures the production 
of a non-empty confidence region. There remains the question 
of what is represented by a possible contention inferable 
from equation (4-5), namely that positive probability p is 
A ached to a single point, the least squares solution. Any 
such statement would be puzzling when approached with the 
viewpoint that specific random variables on a continuum have 
zero probability and that probability statements can only be 
made about boundaries for continuous random variables s 
as location coordinates. 

An answer may be attempted by first pointing out that 
azimuth errors are unknown and independent random variables, 


but they are attached to actual and unknown emitter azimuths 
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with respect to the set of intercepting DF sites. Observed 
femal azimuths have two components, the first of which, 
emitter azimuth, is dependent upon emitter location. The 
second component, azimuth error, is independent from site to 
Site and is the theoretical basis for the Chi-Square contour. 
EEUdence region. Prior to the production of a particular 
confidence region, the normality assumption permits the 
Specification of upper bounds for the Chi-Square statistic 
corresponding to various levels of confidence. During 
peesuction of the particular confidence region, the require- 
ment that a single location must satisfy the observations 
places a lower bound on the Chi-Square statistic. 

Each emitter azimuth can depart from the observed signal 
azimuth, but there are conceptual bands (based on the error 
distribution assumptions) bounding such departures for any 
given probability. When more than two such bands are over- 
lapped in the solution region, then a proportion of the 
possible combinations of such departures are ruled out by 
the constraint that azimuth errors are attached to emitter 
azimuths. Whenever f* is greater than zero, then the prob- 
ability, denoted p*, associated with f* as an upper boundary 
from equation (4-5) represents the proportion of impossible 
sums of squares out of those permitted by the density func- 
tions governing azimuth errors. 

A thought experiment may clarify these remarks. Suppose 
it were possible to specify all combinations of azimuth 


errors permitted by their respective density functions prior 
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to finding the least squares solution. There exists a loca- 
tion satisfying one of the postulated combinations of azimuth 
Ennors for a specific set of signal azimuths only if applica- 
tion of each of the postulated azimuth errors in the combina- 
tion to its corresponding signal azimuth yields a corrected 
set of emitter azimuths whose corresponding great circle arcs 
intersect at a point. Few such postulated combinations in 
the infinite set can be satisfied since the point of inter- 
section for any two corrected emitter azimuths uniquely 
determines the required azimuth corrections for the rest of 
the Observations. Furthermore, there cannot exist a satis- 
meng location for any postulated combination of azimuth 
errors having a sum of squares less than f* since no such 
Corrections are of sufficient magnitude to result in a 

common intersection. 

These suggestions can be illustrated by an extreme and 
improbable example of emptiness of the confidence region for 
emee signal azimuth great circle arcs forming a cocked hat. 
Suppose that the relationships between the DF sites and the 
emitter are as given in Figure 8, and suppose that the errors 
associated with the three DF sites are all governed by the 
standard normal distribution (variance equal to one degree 
squared). Suppose also that each of the three signal 
azimuths were in error by minus two degrees (where the 
probability that azimuth error is less than or equal to 
minus two degrees is equal to .0227) in order to form one 
of the two largest possible cocked hats for о; | equal to 


two degrees. Then the least squares solution will be found 
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FIGURE 8 
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mî tie Center Of the cocked hat, and f* is about equal to 
twelve degrees squared. For three degrees of freedom, there 
E0992 probability that the Chi-Square statistic would 
be less than or equal to twelve degrees squared, indicating 
an empty confidence region at the .992 level (or any lower 
levels). Now the probability that |o. | 15 greater than four 
degrees is equal to .000, suggesting that a boundary of four 
degrees on either side of each signal azimuth great circle 
arc should contain all possible emitter azimuths from which 
each signal azimuth represents a departure. The situation is 
imebustrated in Figure 8. 

With each emitter azimuth free to be anywhere within its 
dashed band by the normality assumption, but dependent upon 
the actual emitter location and therefore constrained by the 
other dashed bands, the shaded intersection of the three 
bands must (for any non-negligible probability) contain the 
emitter. Moreover, the Chi-Square contour methodoloay would 
produce a circular confidence region vithin the shaded 
triangle for any confidence level greater than .992. The 
center.of the shaded region corresponds to the least squares 
solution. Suppose a "center" to be the least squares solu- 
tion united with all points forming a neighborhood within 
radius e of the center. Then this "center" had a probability 
s1 9392 F f(e) of containing the emitter location. Thus 
there was a probability of .008 - f(e) of finding the emitter 
within a region bounded by a "centerless" circle inscribed 


in the shaded triangle. There is zero probability that the 
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emitter lies outside such an inscribed circle, or confidence 
region at the 100 percent level. The description of the 
Fon tclidence region as circular applies to Figure 8 only 
because the DF sites are equidistant and symmetric around 
Ege emitter, and because azimuth errors are minus two degrees. 
ane implication of the foregoing is that if errors are 
normally distributed with known variance, then abnormally 
large cocked hats would be indicative of excellent solutions 
via the least squares approach. The point would be difficult 
for operators to accept since high confidence is normally 
associated with a set of great circle arcs corresponding to 
Signal azimuths that nearly intersect at a point, and large 
departures are suspect. The discrepancy between the above 
conclusion and experience may be attributed to the fact that 
abnormally large cocked hats only occur when outliers (signal 
azimuth errors that do not obey the empirically estimated 
distribution) enter the data base. 

Sometimes strange relationships between the sizes of 
cocked hats and accuracy of locations (least squares solu- 
tions) can also be demonstrated with a network of three DF 
Sites with errors governed by a Bernoulli distribution, 
having equal probability that azimuth error is plus or minus 
three degrees. Figure 9 illustrates the situation. Given 
EN utr ubtoshescenter of the Star, observed signal azi- 
muths are in error by plus or minus three degrees, each 
mith probability 1/2. All possible configurations of 


cocked hats are therefore represented on the figure for the 
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sketched relationship between DF sites and emitter location. 
Six of the possible cocked hats are the small triangles at 
jm@emtips of the star. The other two possibilities are the 
union of the body of the star with the small triangles at 
Pores NE, SE and W, and the union of the body of the star 
with the small triangles at points E, SW and NW. The latter 
two cocked hats are abnormally large and have £* equal to 27 
degrees squared (which is over twice the value of a Chi-Square 
EuUUStic at the .995 level for a normal error distribution). 
Yet, the center of these largest cocked hats coincides 
approximately with the emitter location, indicating excellent 
least squares solutions, and the centers of the triangles at 
the points of the star, the smaller cocked hats, are far 

from the emitter location, mU though their f* value is 
relatively small. 

The probability of occurrence of such large cocked hats 
when error is governed by a normal rather than Bernoulli 
distribution is negligible, and they do not appear when data 
is simulated with a normal error distribution having reason- 
able variance. However, in about 400 simulations during 
development of the Fix Program they were forced to appear 
by governing error with such a Bernoulli distribution. As 
would be expected from Figure 9, a plot of location error 
versus f* verified that location error decreases as f* 
approaches its maximum of 27 degrees squared for this unusual 
Peron distribution. Such a simulation naturally yields about 


25 percent of the abnormally large cocked hats and 75 per- 


46 





cent triangles at the points of the star. With a random 
distribution of DF sites and signal emitters over the globe, 
the stars are not so well-shaped as that in Figure 9. 

If one regards the fix problem with least squares solu- 
d confidence region construction at confidence level 
a as a Bernoulli trial with success if the confidence region 
mains the emitter location and failure if the emitter is 
outside the confidence region, then the expected number of 
successes for m independent trials is am. That the Chi- 
Square contour confidence regions have normalized frequency 
Containment about equal to the probability of containment 
is therefore expected to be true and was demonstrated to be 
true by simulation. At confidence level a, on about a per- 
cent of the trials non-empty confidence regions are produced 
which contain the actual emitter. On the remaining (1 - a) 
percent of the trials, outcomes are distributed one empty 
regions and non-empty regions with the actual emitter out- 
side the confidence region. The lower the confidence level 


the higher the proportion of empty regions amongthe failures. 


D. FORTRAN CONFIDENCE REGION PROGRAM 

A plot subroutine is included in the Fix Program of this 
paper. Sample plots are included in the Computer Output 
section of the paper. The plot is a square grid of constant 
dimension, scaled in nautical miles to contain the confidence 
region. An asterisk on the grid indicates that the corres- 
ponding point has an objective function value greater than 


the Chi-Square boundary value for the specified confidence 
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level. The least squares solution is at the center of the 
m ала iS indicated by the letter S. During program 
development, it was found that a 90 percent confidence level 
is sufficient to preclude the occurrence of empty confidence 
regions, and the sample regions are of that size. 

| ШОО Ртодгаш Description 

Given in a data statement of the subroutine аге 
Chi-Square values corresponding to the numbers of possible 
degrees of freedom, up to n for a network of n DF sites, and 
the selected confidence level. The confidence region boun- 
dary value is selected according to the number of signal 
azimuths available to find a solution. 

The dimension of the grid is 41 rows and 51 columns, 
tailored to the line printer distances between rows and 
columns, so that the ratio of the distances between rows and 
columns is 4:5. Depending upon the scale of the grid, coor- 
dinates on the globe are calculated to match each of the 
ESI points on the grid. Curvature of the earth calculations 
are used to determine distances between points on the grid, 
thereby avoiding distortion inherent in planar assumptions. 

At each of the points on a row, the sum of squared 
azimuth differences is computed for comparison with C 


(n,a)* 
an asterisk being assigned for grid entry if Ir Е) 15 

f 
exceeded. After the confidence region is exited on a given 
row, remaining points on the row are assigned asterisks 


without computation. Similarly, after the confidence region 
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is completed, the rest of the grid is filled with asterisks 
EuSEhout computation. 

The area of the confidence region is estimated by 
counting the blanks and multiplying their sum by the area 
represented by each cell, 80 square miles in the case of 
a 400 by 400 nautical mile grid. If computations on the 
first row yield any points with objective function value less 
than or equal to EC then the grid is rescaled to cover 
a greater area, and the computation starts over. If the 
confidence region leaks out on either of the sides of the 
grid during computation, then again the grid is rescaled, 
EEComputation starts over. In practice, the initial area 
covered by the grid can be tailored to the network and 
expected size of confidence regions. 

2. Results 

ШЕКЕ е Solution search in the Fix Program should 
terminate in the region of the least squares solution before 
the least squares point is found, then the confidence region 
is noticeably asymmetric about the center at S. The lack of 
symmetry indicates that the fix solution should be recomputed 
with a smaller termination step. The center of the confi- 
dence region, or least squares location, could also be 
estimated from the grid itself, but only to an accuracy 
limited by the dimensions of a cell. 

The confidence region subroutine was used successfully 
in over a thousand runs with simulated data. Its chief 


deficiency is the computer time required to produce the grid 
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entries, averaging slightly over ten seconds for numbers of 
signal azimuths per target varying uniformly from three to 
nine. On a slower machine than the IBM 360/67, the subrou- 
tine would require revision to be usable for production. 
Sampling every other or every third cell would be possible 
avenues for reducing computation time, but the present 
version was preferred for a more accurate estimate of con- 
fidence region area. A more sophisticated sns гог reduc- 
ing computation time would be to use a modified Golden 
Section Search to find the confidence region boundary outward 
from the center of the region, store coordinate pairs of the 
boundary in an m by n array, and then print the boundary on 


a scaled grid. 
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ENDE RACER GEOMETRY 


The Fix Program produces several statistics that are 
useful in describing intercept geometry. The first group 
Maese statistics is related to signal azimuths and is 
used to plot position lines as linear approximations of the 
Meat circle arcs in the region of the solution. In the 
second group are some values to which location error is 


memcitive. 


A. PLOTTING POSITION LINES IN SOLUTION REGION 

Following the optional summary of search progress, the 
least squares solution and the value of f*, the computer 
output has a table of values in which each row corresponds 
to one of the intercepting DF sites. Row entries include 
the site identification number, the site name, the azimuth 
Acne solution, the signal azimuth, their difference, the 
range in nautical miles from the solution to the site, and 
three values used > plot the position line relative to the 
least squares solution. These values are calculated by 
ПЕ Лоу with a right spherical triangle illustrated in 
Figure 10. 

Figure 10 is labelled to correspond to the base of the 
polar triangle of Figure 3 with the DF site at vertex B and 
the least squares solution at Vertex C. Vertex angle A is 
tant angle formed by connecting C to the closest point of 


approach (CPA) of the great circle arc related to the signal 
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azimuth 0. Arc a, connecting the DF site and least squares 


solution, is related to the solution azimuth, denoted f*, 
and is given by equation (2-11), using the polar triangle of 


Figure 3. By the law of sines for spherical triangles, 
sin b = sinúu Sin a (5-1) 


The curvature of earth distance covered by arc D is 
EHoOSsely approximated by multiplying arc b in radians by the 
mean earth radius of 3439.9 nautical miles, yielding the 
nautical mile departure of the CPA from the solution. 

The site azimuth with respect to the solution, denoted 
Es found after calculating angle C of Figure 3 in the 
Same manner that vertex angle B and azimuth B were calculated 
in equations (2-12) and (2-13). Then the azimuth of the CPA 
with respect to the solution is the sum of Y* and C, in the 
case illustrated. Finally, the azimuth of the DF site with 
respect to the CPA is found by subtracting 90 degrees from 
the sum of y*and C, in the case illustrated. Similar argu- 
ments apply to three other general cases, solution east or 
west of the DF site and w positive or negative, and resulting 


conditional equations are included in the Fix Program. 
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The azimuth and distance of the CPA with respect to the 
solution are the components of a vector which can be plotted 
on polar coordinate paper. A line of position perpendicular 
to this vector approximates the great circle arc in the 
region of the solution. The direction of the intercepting 
DF site is given by the azimuth value in the last column of 
the output row. The plot can also be constructed with a 
parallel ruler by orienting the ruler along the site azimuth 
with respect to the CPA, moving the free edge the distance 
Seale OL the CPA toward the azimuth of the CPA. Such 
lines of position are plotted to the grid scale on several 
of the confidence region plots provided in the computer out- 


Sut Section of this paper. 


КО MEAN ACUTE ANGLE OF INTERSECTION 

There are factors that may be said to reduce location 
error, given that signal azimuth error is present but un- 
known. One such factor is the number of DF sites providing 
Signal azimuths to make a fix: the more independent measure- 
ments, the TT the emitter location estimate by least 
squares Solution. Another factor is the layout of the DF 
network, the obvious supposition being that a network spread 
over a large area produces better locations of globally 
distributed targets than a network confined to a small area. 

A slightly less obvious statistic is the magnitude of 
the acute angle of intersection of great circle arcs corres- 
ponding to signal azimuths. For two DF sites participating, 


the dependence of location error on the acute angle of 
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intersection is illustrated in Figure ll. In case a, two 
lines of position intersect at about 90 degrees; in case b, 
the acute angle of intersection is about ten degrees. The 
difference in the sizes of regions contained by the dashed 
lines demonstrates the effect of two degree uncertainty of 
measurement on either side of a signal azimuth on the two 
different intersections. 

Although the confidence regions described in Section IV 
of this paper do convey this dependence, that for case a 
being a near circle contained in the dashed lines and that 
for case b being an elongated ellipse of orientation and 
dimension similar to to the dashed region, the acute angle 
An tersection statistic was found to be significant enough 
to warrant isolated analysis. It can be demonstrated 
graphically that the dependence also exists for n greater 
than two. For all cases, the dependence Can be examined by 
calculating the mean value of all the acute angles of inter- 
section formed by all possible permutations of two lines of 
position, given by (n)>- The mean acute angle of inter- 
section is part of the Fix Program output. Dependence of 
location error on this statistic is examined in later sec- 
tions of this paper. For individual pairs, the acute angle 
of intersection is closely approximated by taking differences 


from the site azimuths with respect to the CPA. 
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Sen, MEAN ARC STATISTIC 

enore Subtle factor affecting location error is a 
function of arc length or distance from the DF site to the 
emitter location. The arc from the DF site to the emitter 
and all other arcs representing possible departures due to 
erroneous Signal azimuth measurements converge at the DF 
Site and at its antipode. Therefore, the maximum nautical 
mile departure of a line of position from the emitter loca- 
ESM гот constant azimuth error occurs at 90 degrees of arc, 
mS bout 5403 nautical miles, from the DF site. 

ime arc from the DF site to the least squares solution 
1S given by equation (2-11). The arc statistic is equal to 
this arc, or to its supplement if the arc is greater than 
90 degrees. The dependence of location error on the arc 
statistic is illustrated in Figure 12, where the solid line 
represents an arc from the DF site through the emitter loca- 
tion to the DF site antipode, and the dashed lines represent 
three degree departures from this arc. 

The distance from the site-to-emitter arc to the depar- 


ture arcs increases non-linearly from zero at the DF site 


c ___ _>__ _________J__J_____ = ===—— 
Site $ Antipode 


Relating Error to Emitter Distance 


FIGURE 12 
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Ets antipode to 180 nautical miles at the midway point. 
EN later section of this paper, location error sensitivity 
ES the mean arc statistic for all DF sites participating in 


mem location is examined. 


S 








DIS IMUDIDIONSPROCGEDUREE 


A modified version of the Fix Program was used to simu- 
late operation against targets. Targets were generated by 
a subroutine which provided longitude and geocentric lati- 
BEnuesdescribing random locations uniformly distributed over 
the globe. Such coordinates vere saved for later comparison 
with the least squares solution. 

The number of intercepting sites for each target was 
either specified for simulations on fixed degrees of freedom, 
Sr Tandomly generated according to a uniform distribution 
from three to nine. In both cases, the specific intercept- 
ing DF sites were randomly selected from a set of 19 DF sites 
composing the particular network being simulated.  Azimuths 
pene target location with respect to the participating DF 
sites were then calculated. For each such emitter azimuth, 
a random error was added, with direction and magnitude of 
the error governed by the standard normal distribution. The 
emitter azimuth with error was then treated as a measured 
Smal azimuth for input with the site identification number 
to the Fix Program. | 

During program development, the simulation procedure was 
used with zero signal azimuth error generated to insure that 
the least squares solution of the search procedure coincided 
with the target coordinates to at least an accuracy of three 


decimal places. A search termination value of one-tenth of 
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a nautical mile was found to be adequate, and the present 
combination of Golden Section Searches and second-order 
gradient searches was found to be satisfactory during the 
development simulations. The Fix Program has also been 
ENS esstully tested with up to 19 participating DF sites on 
a single target, with a set of DF sites contained in a 180 
by 240 nautical mile area, and with higher variance error 


miserioutions, including outliers. 


A. GEOCENTRIC AND GEOGRAPHICAL LATITUDE 

A comparison of least squares solutions using geocentric 
and geographical latitude in Fix Program calculations was 
made in order to examine the accuracy gain of the geocentric 
model. Errorless signal azimuths of random targets were 
computed using geocentric latitude in calculations. A least 
squares solution was then found by using geographical lati- 
tude in calculations, and the conversion of the solution to 
geographical latitude was omitted. This distance between 
the geographical coordinates of the target and the pseudo- 
geographical coordinates of the least squares solution 
measures the departure from the least squares location due 
to use of geographical latitude on an assumed sphere. 

Radial departures for two degrees of freedom were occa- 
sionally devastating (424 nautical miles in one case). Using 
an imaginary network, called the national network, and degrees 
of freedom uniformly distributed from three to nine, the mean 
radial departure for 148 samples was 6.8 nautical miles. 

The direction of such departures tends toward the equator. for 


targets in both northern and southern hemispheres. 
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NANI LOCATION ERROR ANALYSIS WITH FIXED DEGREES OF FREEDOM 


The sensitivity of location error to Several factors was 
examined in simulations with degrees of freedom fixed from 
two to seven. From 100 to 200 samples were generated at each 
level for each of the networks. All of the following plots 
are drawn to the same scale to simplify between plot compari- 
sons of the effects of different degrees of freedom and 
different network configurations. Values too large for the 
common scale are entered in a row at the upper limit of the 


plot above an off of scale indicator. 


A. COMPARISON OF ERROR ON NATIONAL AND GLOBAL NETWORKS 

Member sites of the national network varied in longitude 
From Attu to Boston and in latitude from Anchorage to Hono- 
lulu. The specific sites are identified in the second data 
Statement of the Fix Program. Sites in the widely distri- 
buted global network included Anchorage, Bombay, Boston, 
Buenos Aires, Cape Town, Christchurch, Dakar, Honolulu, Houston, 
Lima, Naples, Papeete, Perth, Recife, Seattle, Seychelles, 
Tehran, о ања Wake. 

Figure 13 is a plot of mean location error from each of 
the sample groups against the number of intercepting sites. 
Dashed lines connect the means of each network. The plot 
shows nearly double the location error in the national net- 


work even at higher degrees of freedom. 
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With all other factors random, the average accuracy of 
a network is dependent upon the network layout and the 
variance of signal azimuth measurement error. For a given 
network, the mean location error would vary in the same 
direction as changes in such variance, in all of the follow- 


ing simulations held constant at one degree. 


B2 ERROR VERSUS MEAN ACUTE ANGLE OF INTERSECTION 

Based on samples from the national network, Figures 14, 
land 16 are scatter plots of location error versus the 
mean acute angle of intersection explained in section V-B. 
Degrees of freedom were fixed at three, five and seven, 
respectively. 

The strong correlation in the national network can be 
interpreted as a reduction of both the mean and variance of 
location error as intersections of great circle arcs are less 
Beute. Rough scatter plots for two, four and six degrees of 
freedom showed the same correlation, naturally most pro- 
nounced at two degrees of freedom. The set of plots shows 
a shrinking of И Жанес of the mean acute angle of inter- 
section statistics as degrees of freedom are increased. 

Figures 17, 18 and 19 are similar scatter plots with 
entries from simulations on the global network. The corre- 
lation is still evident at three degrees of freedom but be- 
comes masked at higher levels as the layout of the mm 
reduces the chance of producing a combination of target and 


intercepting sites with a small mean acute angle of inter- 
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section. Consequently, shrinking of the range of the 


statistic is extremely evident on Figures 18 and 19. 


ERROR VERSUS MEAN ARC STATISTIC 

Figures 20, 21 and 22 are scatter plots of national 
network location error versus the mean arc statistic 
(explained in section V.C) for three, five and seven degrees 
of freedom, respectively. At all levels, the large location 
errors tend to be associated with the high end of the range 
of mean arc statistics. The strength of the correlation 
appears to increase as degrees of freedom increase on 
succesSive plots. 

Figures 23, 24 and 25 are similar plots for the global 
network. Again, the correlation is masked after three 
degrees of freedom with the shrinkage effect extremely pro- 


nounced as the degrees of freedom increase. 
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ED c LOCATION ERROR ANALYSIS WITH VARYING DEGREES OF FREEDOM 


Simulations were conducted on both networks with degrees 
Gf freedom generated uniformly from three to nine. Chi- 
square contour confidence regions were generated for each 
target's least squares solution. Using a confidence level 
of 90 percent, national network confidence regions contained 
90.38 percent of the generated targets. Results for the 
global network were 90.44 percent. Such close agreement 
between the confidence level and containment in the simula- 
tions is dependent upon a reliable standard normal azimuth 
ET generator. In production use, such reliability could 
only be achieved by having a good estimate of DF site 


Variances. 


A. NATIONAL NETWORK 

Figure 26 is a scatter plot of confidence region area 
versus location error. Again using a common scale to aid 
network MEE seven points off the scale are entered 
above or to the right of an indicator. Although there is 
one hopeful cluster near the plot origin, the large variance 
location errors at given confidence region areas is 
maerecative of the futility of predicting location error 
based on the size of this reliable confidence region model. 

Kalte 27 is a multi-level histogram of location error 
at the various degrees of freedom, demonstrating the same 


Bout Of correlation evident in Figure 13. Figure 28 
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een trates the correlation between location error and the 
mean acute angle of intersection, present even with mixed 
degrees of freedom. The same holds true in Figure 29, a 


fot Of location error versus the mean arc statistic. 


B. GLOBAL NETWORK 

Figures 30 through 33 are similar plots for the global 
network. On Figure 30, location error and confidence region 
area appear independent. Figure 31 shows good correlation 
between location error and degrees of freedom. Figures 32 
and 33 show nearly complete masking of the previously 


established correlations due to the layout of the network. 
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О CONCLUSION 


mac Fix Program and Chi-Square contour confidence region 
Subroutine could be integrated into the software of position 
location systems by direction finding. The present program 
meen target and azimuth error generators included) could be 
used to simulate an existing DF network for accuracy studies 
on potential targets in different regions of the globe or for 
studies on optimal placement of DF sites in a proposed net- 
work. 

The relationships demonstrated in the sections on loca- 
tion error analysis indicate that larger location errors 
accompany contraction of the network layout; however, this 
factor would be mitigated by a practice of assigning regional 
rather than globally distributed targets to the network for 
prosecution. 

Countering the widespread belief that small confidence 
regions indicate accurate locations and large confidence 
regions indicate poor locations, it has been demonstrated 
that for a reliable model location error and confidence 
region area are nearly independent. 

The relationships examined as parts of location error 
analysis suggest an alternative answer for the question of 
how accurate is a given location. For an existing network, 
Mata could be collected from locations of known transmitters. 


Such data would include statistics on number of signal 
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azimuths available to make the fix, mean acute angle of inter- 
section and mean arc statistic. With known location error 

as the dependent variable, multiple linear regression 

analysis could be used to develop regression and correlation 
coefficients for predicting location error as a function of 
independent variables measured or calculated at the time of 

lo ition. 

The workability of such an approach was not demonstrated 
in the present paper for several reasons: I The effort 
required to perform the multiple linear regression (using a 
standard program from the IBM scientific subroutine package) 
would be a project in itself and would be best reserved for 
actual data from an existing network. (2) System operators 
would be aware of other measurable factors for inclusion as 
independent variables in the regression analysis, such as 
variance on given signal azimuth measurements, signal 


Strength and measurement duration. 
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